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Abstract 

Three ways of fitting the tt + tt~ mass spectrum in y(4260) decays are studied. Data 
presented recently by Belle cannot be fitted by the tttt S-wave intensity for elastic scattering. 
They can be fitted by adding a rather arbitrary destructive interference with the a pole term. 
A better fit may be obtained with the decay sequence Y — > ttH, H — > irJ/ip, where H is 
a J PC = I— ccqq state peaking at 4.0 GeV with a width ~ 280 MeV. A third possibility, 
involving a triangle diagram due to y(4260) -> Z)L>i(2420), D 1 -> Z?*(2007)vr, DD* -> vrJ/^ 
fails to fit the data. The first and second possibilities (or a combination) could be resolved 
by analysis of the Y(4260) Dalitz plot and nJ/ifj mass projection, which are not presently 
publicly available. 

PACS: 13.25.-k, 13.75.Lb, 14.40.Cs, 14.40.Ev. 
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There has been considerable interest in the V(4260) discovered by BaBar [1]. It has been 
5 confirmed by Cleo [2], by further observations by BaBar [3,4] and Cleo [5] and most recently by 
^ .Belle [6]. Its clear observation in decays to TnrJ/ip has led to suggestions that it may be related 
to the ccg hybrids predicted in this general mass range [7-10]. It is observed in the radiative 
return process, where it is produced through e + e~ — > 7; it therefore has J PC = 1 . Cleo observe 
V(4260) — > 7t°7t°J/iIj [2] with a rate consistent with half that of 7T + 7i~J/ip, albeit with just 8 clean 
events. This result strongly supports / = for F(4260). 

The latest Belle data [6] reveal a curious ttit mass spectrum. The objective here is to examine 
• 1 ways of fitting this spectrum. Three classes of decay scheme for F(4260) are investigated. 

The simplest hypothesis is Y — > (7rn)sJ/ip, where (vrvr)^ stands for the S-wave as observed in 
o3 elastic scattering. The full curve of Fig. 1(a) compares data points with the spectrum predicted 
from the latest parametrisation [11]; a brief subroutine for this parametrisation is available from 
the author. In making this prediction, the line-shape of F(4260) has been folded with the phase 
space for Y — > (7111)3 J /ip using the Breit-Wigner parameters determined by Belle. 

The predicted peak at ~ 900 MeV is too wide, fails to reproduce the low point at 750 MeV 
and does not account for the low mass peak. A minor detail is that none of the fits reported here 
reproduce accurately the highest mass point at M{jnx) ~ 1275 MeV. This mass is outside the 
range accessible for Belle's central Y mass, 4295 MeV; it is possible that events in this bin arise 
from spreading of the high mass tttt peak due to mass resolution. This bin is therefore omitted 
from all fits, and the curve is normalised to the 113 events in the remaining 9 bins. The resulting 
fit has x 2 = 69 for 9 bins and 8 degrees of freedom. 
In elastic scattering, the tttt amplitude is 

f(mr s N ( s ) MT ^) m 
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Figure 1: (a) The hit mass spectrum from Belle data for F(4260) — > 7r + 7r~ J /ip (data points), com- 
pared with (i) production of J/ip with the elastic ttit S-wave (full curve), (ii) including destructive 
interference with the a pole using eqn. (2) (dashed curve), (iii) with the triangle diagram of Fig. 
3 (dotted); (b) the corresponding comparison with Cleo data. 



where T(s) is close to linear in s, with a zero at the Adler point s = sa = [The convention 

adopted here is that the J ftp is taken as particle 3, hence s i2 is the invariant mass of nn]. The 
-kit phase shift goes through 90° at M~ 950 MeV and the intensity of irir elastic scattering peaks 
there. However, the pole in the amplitude lies much lower, at ~ 470 MeV, because of the strong 
s-dependence of the width. In some production processes, e.g. J/ip — > unix [12] and D — > 37r 
[13], the production amplitude is instead well fitted by the same Breit-Wigner denominator D(s) 
as eqn. (1), but with a constant numerator N(s). The elastic amplitude seems to give a good fit 
to processes involving low momentum transfers, e.g. decays ip(2S) — > irnJ/ip [14]and — > 77r°7r° 
[11]. The production processes of Refs. [12] and [13] involve larger momentum transfers, hence 
more distant singularities. The decay V(4260) — > TrnJ/ip is closer to the former case. 

An adequate fit may be obtained using destructive interference between Eqn. (1) and the a 
pole by altering the numerator N(s) to 

, (g - sa) - G(M 2 - sa) 

fs(**) = ^ , (2) 

with the same Breit-Wigner denominator D(s) as in eqn. (1). The resulting fit is shown by the 
dashed curve in Fig. 1, with G = 0.37, x 2 = 17.8 for 9 bins (7 fitted parameters). This fit seems 
somewhat artificial, since the numerator has no obvious explanation; however, it cannot presently 
be ruled out. Making G complex does not improve the fit significantly. 

It seems unlikely that / (980) is responsible for the peak at 950 MeV in Fig. 1, because the f 
is too narrow. The latest BES data on J/ip — > 07T7T and <pKK [15] determine the full- width at half 
maximum of /o(980) — > tttt to be 34 ± 8 MeV. Belle do not quote their 7ttt mass resolution, but 
do quote the width of V(4260) as 133 ± 26^ 3 MeV. The mass resolution in both cases should be 
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Figure 2: (a) Fit to the Belle hit mass spectrum for L = 1 decays via an intermediate H with 
M = 4080 MeV, T = 280 MeV; (b) the upper half of the Dalitz plot; (c) fit to Cleo data. 



comparable. In any case, the 7nr mass distribution cannot be fitted by /o(980) and the a pole: the 
contribution from the a pole is too broad if it interferes with only the very narrow / (980). Cleo 
observe 2 events for the decay to K + K~ J/ip [2]. These could arise from the known amplitude for 
a — > KK or from a small / (980) J ftp contribution. 

Fig. 1(b) compares the same fits with the 37 events of Cleo [2]. Although statistics are lower, 
the peak in data again looks narrower than the full curve; there is also some tentative indication 
of a rise at low tttt mass. Cleo correct for background from sidebins around the J/ip, but do not 
correct for possible background under the V(4260) from other if) states. Judging from the mass 
range 4.0 to 4.2 GeV, this background could be 15%. 

BaBar [1] show a nn mass spectrum which does not show any rise at low ixix mass. It has 
been obtained by summing the V(4260) signal over the whole mass range and then making a 
background subtraction. The resulting signal/background ratio is ~ 1 : 1. In this respect, Belle's 
strategy of selecting events in the mass range 4.2 to 4.4 GeV, i.e. 1.5 times the full width, gives 
a better signal/background ratio for the purpose of viewing the nn mass spectrum. Some such 
weighting of data points near the peak of V (4260) is desirable. 

Let us turn to a second hypothesis: sequential decays of the form Y — > ttH, H — > nJ/ip. If Y 
has 1 = 0, as Cleo data suggest, the H has 1 = 1 and could be a 4-quark state. Two cases need 
consideration. When the orbital angular momentum L in both steps of the sequence is zero, H 
has J PC = l + ~. UL = 1,H has J PC = 1", 0"" or 2—. 

With L = in both steps, the predicted nir mass spectrum is close to the full curve of Fig. 
1 and hence fails to fit the data. However, with L — 1, a good fit emerges naturally, as shown 
on Fig. 2(a). This result depends on simple kinematics. Consider the case V123 — > 
H23 — > ^2{J where numbers label the particles. If L = 1 in the second step, there is a factor 
cos 62 in the amplitude; here 62 is the angle of 7r 2 in the 23 rest frame with respect to the recoil 
direction of particle 1. For fixed s 2 3, s 12 varies linearly with cos 6*2 from one side of the Dalitz plot 
to the other, leading to a zero halfway between left- and right-hand edges of the plot. The mass 
of H determines the value of S12 = m\ 2 a ^ which the zero occurs, via the boundaries of the Dalitz 
plot. The width of H controls the sharpness of the dip in the mass spectrum near 750 MeV. 
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Fig. 2(b) shows the upper half of the Dalitz plot for this hypothesis. The zero due to the 
factor cos 9 2 is clearly visible. Fig. 2(c) shows the fit to Cleo data. The fitted mass and width 
are almost linearly correlated because of the L = 1 centrifugal barrier. The optimum fit has a 
mass of 4.08 GeV with T = 0.275 GeV and x 2 = 15.1 for 2 fitted parameters, seven degrees of 
freedom. Fits poorer in x 2 by 1 are obtained with M = 4.04 GeV, T = 0.225 GeV or M = 4.16 
GeV, T = 0.330 GeV. Because of the centrifugal barrier for production, the peak intensity is just 
below 4 GeV for all solutions. A Y mass at the lower end of the range is most plausible, in order 
to give a reasonable separation in masses of Y and H. 

In detail, the amplitude for production of H 23 alone is 

/ oc /y(si 23 )/i?(s23)-Bi(Pi)-Bi(P2) cos 6» 2 , (3) 

where /y is the Breit-Wigner amplitude of Y as a function of S123, fu is the Breit-Wigner amplitude 
for H, Bi(pi) is the L = 1 centrifugal barrier factor as a function of the momentum of 7Ti in the 
rest frame of Y, and Bi(p 2 ) is the centrifugal barrier factor as a function of the momentum of n 2 

in the rest frame of H 23 . The centrifugal barrier is taken as Bi(p) = pj \jp 2 + p^ where p = 0.25 
GeV/c. The centrigual barrier cuts off the intensity near the top of Fig. 2(b). The amplitude of 
eqn. (3) needs to be symmetrised between particles 1 and 2. The usual weighting according to 
phase space ds\ 2 ds 23 is also applied. 

None of the experimental groups has presented data as a function of m 2 (n J / tp) . Even with the 
current low statistics, structure like that of Fig. 2(b) should be clearly visible. It gives rise to a 
strong peaking at two corners of the Dalitz plot and a peak in the nJ/ip spectrum near 4.0 GeV. 
The situation could however be confused by interferences between amplitudes of the two forms 
considered so far. Furthermore, interference between V(4260) and underlying background from 
^(4160) and -^(4415) is possible. Indeed, van Beveren and Rupp have argued against the existence 
of F(4260) as a separate resonance [16]; they interpret the structure in this mass range as arising 
from thresholds and interferences in the s-dependence of -^(4040) and -^(4160) amplitudes. 
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Figure 3: The triangle diagram for Y -> DZ^^O), D x (2420) -> ttD*(2007), DD* (2007) -»• 

iyJ/i/j. 

A third hypothesis will now be considered. The V(4260) lies close to the threshold for 
-D(1865)-Di(2420), particularly for the resonance mass quoted by Belle, 4295 MeV. A possible 
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mechanism for production of the imJ/ifj system is shown in Fig. 3. In this so-called triangle dia- 
gram, the Di decays (via a D-wave) to 7rD*(2007). This long-lived resonance can de-excite from 
the spectator D producing nJ/ip. A simulation of this process gives the dotted curve of Fig. 1(a), 
and fails to fit the Belle data. Even if interference with the (n7r)s J/i/j amplitude is included with a 
fitted complex coupling constant, nothing resembling the Belle mass spectrum can be reproduced. 
The essential feature of the process in Fig. 3 is that the D and D\ are produced almost at rest in 
the Y rest frame. The pion from _Di(2420) decay has a momentum close to 360 MeV/c in the D 
rest frame and hence in the Y rest frame. This makes a clear signature for this process. 

In summary, the mr mass spectrum observed by Belle in F(4260) decays is inconsistent with 
that expected for a final state {im) s J/ip unless the numerator of the amplitude gives rise to 
destructive interference. The simplest, but speculative, hypothesis is the sequential decay to J/t/j 
via an intermediate J PC = 1 4-quark state. This possibility fits the data unaided; on the other 
hand, there is no obvious reason why the F(4260) alone should prefer decays to imJ/ip. It should 
be simple to verify or reject this hypothesis in its simplest form, although interferences with the 
first mechanism could complicate the situation. 
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